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Abstract. We use deep (~120 ks) XMM-Newton data of the M87 halo to analyze its spatially resolved temperature structure 
and chemical composition. We focus particularly on the regions of enhanced X-ray brightness associated with the inner radio 
lobes, which are known not to be described very well by single-temperature spectral models. Compared to a simple two- 
temperature fit, we obtain a better and more physical description of the spectra using a model that involves a continuous range 
of temperatures in each spatial bin. The range of temperatures of the multiphase gas spans ~0.6-3.2 keV. Such a multiphase 
structure is only possible if thermal conduction is suppressed by magnetic fields. In the multi-temperature regions, we find a 
correlation between the amount of cool gas (with a temperature below that of the surrounding X-ray plasma) and the metallicity, 
and conclude that the cool gas is more metal-rich than the ambient halo. In the frame of the assumed thermal model, we estimate 
the average Fe abundance of the cool gas to ~2.2 solar. Our results thus point toward the key role of the active galactic nucleus 
(AGN) in transporting heavy elements into the intracluster medium by uplifting cool, metal-rich gas from the galaxy. However, 
the abundance ratios of O/Si/S/Fe in and outside the X-ray arms are similar, indicating that the dominant fraction of metals in 
the gas halo was uplifted by AGN outbursts relatively recently compared to the age of M87. Our best estimate for the mass of 
the cool gas is 5 x 10** Mq, which probably stems from a mixture of ICM, stellar mass loss, and Type la supernova products. 
~30-1 10 Myr are required to produce the observed metals in the cool gas. Finally, we put upper limits on possible non-thermal 
X-ray emission from M87 and, combining it with the 90 cm radio maps, we put lower limits of around ~0. 5-1.0 /jG on the 
magnetic field strength. 
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1. Introduction 

The impact of active galactic nuclei (AGN) is currently the 
most promising solution for the cooling flow problem in galaxy 
clusters: e.g., Binney & Tabor (1995); Briiggen & Kaiser 



(|2002ll; [Bohringer et al., (2002); ,Birzan et al.,(2004) ; pmma 
|et al. ( 2004[ i ; [Sijacki et aL ( 2007|l and references therein; for 
a review see Pete rson & Fabian ( 2006| l; |McNamara & Nulsenj 
( |2007) . It is moreover invoked to transport heavy elements into 
the intracluster medium (ICM) ( Bruggen||2002t |Omma et al.| 
2004[|Rebusco et al.|20 06 ; Roed iger et al.|2007F and to explain 
the exponential cut-ofF of the bright end of the galaxy luminos- 
ity function (e.g. |Croton et al.|2006[ [Sijacki et al.|2007 



The AGN-ICM interaction is one of the driving motivations 
for deep observations that allow detailed analysis of the prop- 
erties and spatial structure of the intracluster gas. While for a 
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complete understanding of the ICM physics a combined anal- 
ysis of multiwavelength data is usually necessary, most of the 
information about the ICM is gained from observations in the 
X-ray domain in which this hot plasma emits most of its ra- 
diation. The main targets of such detailed analyses are thus in 
particular nearby X-ray bright objects, where both good spatial 
resolution and high spectral statistics can be achieved. 

This makes the hot gas halo of M87, which is the second 
brightest extragalactic X-ray source in the sky located at the 
center of the nearby Virgo cluster at a distance of only 16 Mpc 
(Tonry et al..2001) , an ideal target for studying the AGN-ICM 
interaction in great detail. M87 is known to host an AGN pow- 
ered by the galaxy's c entral supermassive black hole with a 
mass of 3.2 x 10'' Mq ( [Harms et al. 1994) . The AGN jet and 
unseen counteijet are believed to drive the complex large-scale 
system of lobes observed in the radio domain (e.g. [Owen et al.[ 
[2000) , which clearly interact with the X-ray gas. The evidence 
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for this is seen primarily in the spatial correlation between 
the two "inner" radio lobes extending east (E) and southwest 
(SW) of the core of M87 and regions of increased X-ray sur- 
face brightness, also referred to as the E and SW X-ray arms 



(|Feigelson et al. 


1987; Bohringer et al. 1995[ Belsole et al.| 


2001 Young eta 


.|2002(|Forman et al.|2005 2007 1. 



The impact of the AGN on the X-ray emitting gas in M87 
is however not limited to the surface brightness enhancement. 
Also the spectral properties in the E and SW arms are affected, 
these regions being the only part of the M87 halo which re- 
quire more complex multi-temperature spectral modeling while 
the rest of the gas is well approximated by single-temperature 
models ( |Molendi|[2002) [Matsushita eFaLlpOO^) [Belsole et aT 



2001 



et al. 



A possible explanation for this is given by Churazov 



(2001 



who suggest that the radio lobes rise buoyantly 
through the hot plasma, uplifting cooler gas from the central re- 
gion and mixing it with ambient gas at larger radii. If this model 
is correct, then M87 offers a unique chance to study the effects 
of AGN-induced gas mixing and gas transport and the influ- 
ence of these phenomena on the energy balance in the cooUng 
core and on the metal distribution in the gas halo. 

In this paper we use deep XMM-Newton X-ray observa- 
tions of the M87 halo to characterize and understand the tem- 
perature structure and spatial variation of the multi-temperature 
gas in the X-ray arms, and to investigate its correlation with 
the radio plasma. Moreover, we aim to map the distribution of 
various metals and verify the influence of the AGN and the ra- 
dio lobes in distributing and transporting these heavy elements 
from the galaxy center into the M87 hot gas halo. 

The paper is laid out as follows: in Sect. |2] we present the 
data sets and data reduction techniques employed; in Sect.[3]we 
describe the various multi-temperature spectral models used in 
our analysis; in Sect. |4] we discuss the thermal structure and 
spatial distribution of the multi-temperature phases and com- 
ment on possible suppression of the Spitzer heat conduction; in 
Sect. [5] we present the spatial distribution of the metals in the 
M87 gas halo and the abundance patterns of different metals 
relative to Fe in and outside the radiolobe regions; in Sect. |6] we 
estimate and discuss the mass and metallicity of cool gas up- 
lifted by the AGN; finally, in Sect. [7] we give an upper limit for 
non-thermal emission from the radiolobe regions. Our conclu- 
sions are summarized in Sect. |8] We adopt a redshift for M87 
of z=0. 00436, a luminosity distance of 16 Mpc ( )Tonry et al. 
200 l[ l, and a scale of 4.65 kpc per arcminute. 



2. Observations and data analysis 

M87 was observed with XMM-Newton for 60 kiloseconds (ks) 
on June 19, 2000 and re-observed for 109 ks on January 10, 
2005. After removing the periods affected by soft-proton flares, 
the net total exposure time is 85 ks for EPIC/pn and 120 ks for 
each EPIC/MOS detector This deep exposure of such a near 
and bright source allows us to analyze the gas halo around M87 
in great detail and with excellent statistics. 

For the analysis, we used the 7.0.0 version of the XMM- 
Newton Science Analysis System (SAS) and employed the 
standard analysis methods as described in e.g. [Watson et al.[ 
( 2001| l. Point sources were found using the source detection 



algorithm implemented in SAS and removed from the observa- 
tion after a visual check to eliminate spurious detections. The 
spectra from each of the two PN observations were corrected 
for out-of-time events separately. We note that our observations 
were affected by gain calibration uncertainties, in some cases 
the spectra - especially for the PN detector - being blue-shifted 
by up to ~1500 km/s. The MOS detectors were also influenced 
by this shift, but the effect was a factor of ~3 less strong than 
for the PN. To account for this, we artificially blue-shifted our 
spectral models. 

As background we used a collection of blank-sky maps 
(Read & Ponman 2003| l scaled according to the correspond- 
ing exposure times for each detector and each observation. The 
accuracy of this scaling for the present dataset is discussed in 
[Simionescu et al.[ ( [2007] ). 

We used an adaptive binning method based on weighted 
Voronoi tessellations (Diehl & Statler 2006), which is a gen- 
eralization of the algorithm presented in Cappellari & Copin 
( [2003) ), to bin the combined counts image to a signal-to-noise 
ratio (SNR) of 200 in the energy range 0.5-2.0 keV, corre- 
sponding to 40000 counts per spatial bin. The binning was 
constrained to follow four different radio contours of the 90 
cm radio image at 10, 3, 1 and 0.3 mJy (each region between 
two successive radio contours was binned separately). The re- 
gion inside the highest chosen radio contour was not fitted in 
order to avoid possible contamination by the AGN, whose flux 
was particularly high in the second observation causing photon 
pileup in the EPIC detectors ( [Werner et al.p OOBa). The spec- 
tra were fitted in the energy range between 0.4-7.0 keV with 
the various models described in the next section. We fitted the 
spectra from all six datasets (two observations with three de- 
tectors each) simultaneously with the same model parameters 
except for a factor to account for variations in the relative over- 
all normalization between the detectors. For this, the PN of the 
new observation was taken as reference. 



3. Spectral models 

For the spectral analysis we use the SPEX spectral fitting pack- 
age (Kaast ra~et al.|1996) . We fit the observed spectra with two 
models: a simple two- temperature MEKAL model in which the 
abundances of the two thermal components are coupled to each 
other, and a differential emission measure (DEM) model with 
a cut-off power-law distribution of emission measures versus 
temperature (wdem). 

The emission measure i' = / n^nidV (where and n\ are 
the electron and ion densities, V is the volume of the emitting 
region) in the wdem model is specified in Eq. ([T]i following 
from iKaastra et al.[ ( |2004l ) : 



dY 
df 



Aj^l/a <' T <r T 

^ mm ^ ^ ^ max? 

elsewhere. 



(1) 



The emission measure distribution has a cut-off at r,nin - 
cTmax- For a — > 00 we obtain a flat emission measure distribu- 
tion. 
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The fraction of the emission measure of the gas lying below 
a chosen temperature Tniaxcooi to the total emission measure is 
given by: 



DEC (deg) 



f maxcool 



maxcool 



dT 



(2) 



By integrating this equation we obtain a direct relation for 
Fniaxcooi/ ^ as a function of a and c: 



^maxcool ( ^maxcool / ^max ) 



1 



(3) 



The wdem model appears to be a good empirical approxi- 
mation for the spectrum in cooling cores of clusters of galaxies 
(e.g. IKaastra et al.|[2004l [Werner et al.1|2006b[ [dePlaa et al.| 
[2006l [ 

In the spatial bins where either the normalization of the 
cool component for the two-temperature fit or the lower cut- 
off for the wdem model were less than 2cr significant, we fixed 
the temperature of the cool component to 1 keV and the lower 
cutoff to O-STmax, respectively. Unless otherwise sta ted we as- 
sume a hydrogen column density of 2.0 x 10^"cm"^ (Lieu et al. 
[1996}. 

Throughout the paper we give the measured abundances 
relative to the proto-solar values given by 'Lodders ' (2003 ) . The 
recent solar abundance determinations by Lodders (2003 ) give 
significantly lower abundances of oxygen, neon and iron than 



those measured by Anders & Grevesse ( 1989 1. 



4. Temperature structure of the M87 halo 

4.1. Two-temperature models 

Two-temperature models are the simplest next step beyond 
the single-temperature model approximation. Based on the 
first XMM-Newton observation of M87, [Belsole et al.| P001), 
Molendi] ( p002j and ^Matsushita et al.| ( |2002| l have already 



shown that the X-ray bright E and SW arms associated with 
the inner radio lobes are significantly better described by two- 
temperature models with a cool temperature component at 
around 1 keV and a hotter component at around 2-2.5 keV. 
With the new data we are able to analyze in much greater spa- 
tial detail the properties of the cool component. 

We find, in agreement with previous work ( |Molendi|2002[ 
Matsushita et al.| |2002\ that the temperature of the cool com- 



ponent remains quite stable between 1 and 1 .2 keV and shows 
little or no trends with radius or with surface brightness (Fig. 
[T]). For example, the correlation coefficient between the tem- 
perature of the cool component and radius in the E arm is only 
0.198. The temperature of the hot phase meanwhile exhibits a 
fairly regular radial gradient (Fig.|2]i. 

For comparison we also performed a similar two- 
temperature model fit using the APEC model in XSPEC. The 
differences in determining the temperature of the cool com- 
ponent can be seen in Fig |3] In general, we find very similar 
trends from both models, however the APEC model seems to 
give overall higher values than the MEKAL model. This effect 
is more pronounced for values of Tcoolapec > 1 -2 ke V than for 
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Fig. 1. Temperature of the cool component as determined from 
the two-temperature MEKAL fit. Only the bins where the nor- 
malization of the cool component was more than 3cr significant 
are shown. 
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Fig. 2. Temperature of the hot component as determined from 
the two-temperature MEKAL fit. The map exhibits a fairly reg- 
ular radial gradient. 
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Fig. 3. Differences in determination of the temperature of the 
cool component with the APEC (fitted in XSPEC) and the 
MEKAL (fitted in SPEX) models. Errorbars are at the Icr level. 
The equality line Ti^mekal - Ti^apec is overplotted. 



TcooLAPEC < 1 -2 keV. We note that the MEKAL model imple- 
mented in SPEX includes a more recent list of transitions than 
the MEKAL implemented in XSPEC (J. Kaastra, private com- 
munication) and that APEC is not implemented in SPEX, thus 
a direct comparison of the two models is not possible. 

While the temperature of the cool phase does not seem to 
vary strongly in the fits, the normalization of the cool gas shows 
spatial variations correlated with the inner radio lobes. In Fig.|4] 
we plot the fraction of the emission measure of the cool phase 
to the total thermal emission, Fcooi/(^hot + ^cooi). overlaid with 
the 90 cm radio contours (O wen et al.|2000 '). As already noted 
by e.g. ,Belsole et al. (2001) , the spatial correlation is more pro- 
nounced in the E arm, while in the SW arm it seems that the 
radio plasma is twisting around the uplifted arm. 

4.2. The differential emission measure model 

We also fitted a model with a continuous temperature distribu- 
tion where the emission measure is a power-law function of the 
temperature, as described in Sect. [3] As the model by Churazov 
|et al. (2001) also indicates, a multiphase structure with a con- 
tinuous distribution of the emission measure as a function of 
temperature is more realistic because, during the rise of the ra- 
dio lobes, gas at several radii gets entrained and mixed (see also 
Briiggen & Kaiser 2002| l. The power-law shape of the emission 
measure as a function of temperature is the simplest model that 
entails a continuous temperature distribution. 

Our main direct result from the wdem fit is shown in Fig.|6] 
where we plot the map of the inverse slope a of the power-law 
emission measure distribution. Higher values in the map indi- 
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Fig. 4. Map of the fraction of the emission measure of the cool 
component to the total emission measure (hot + cool compo- 
nent) for the two-temperature MEKAL fit. 90 cm radio con- 
tours are overplotted in white. 

cate regions where the emission measure distribution is flatter, 
hence where one expects proportionally more cool gas to have 
been mixed together with the ambient gas. To quantify this 
more exactly, we calculate from the fit results the fraction of 
the emission measure of the cool gas to the total emission mea- 
sure according to Eq. |3] We chose Tmaxcooi as 1.5 keV, which 
is roughly the smallest deprojected value of the ambient gas 
computed by Matsushita et al. (2002 1. Fig. |5] shows the spatial 
distribution of the cool gas fraction, revealing very similar spa- 
tial trends as in Fig. |4] where the fraction of cool gas in the 2T 
model was plotted. Thus, the relative contribution of the cool 
gas to the total emission is similar in both models. 

However, the temperature structure of the cool gas can 
be very different, since the same relative contribution can be 
achieved with a shallower slope (large a) and larger low- 
temperature cutoff or with a steeper slope but a smaller cT^^^. 
In Fig.|7] we plot the map of the lower temperature cutoff cr^ax 
in the regions where we find a more than 3cr significant fraction 
of cool gas (in the other regions the gas is not significantly non- 
isothermal, thus it is irrelevant to speak of a lower temperature 
cutoff). We find lower cutoff temperatures which go down to as 
little as 0.6 keV, thus the wdem model allows for the presence 
of gas much below the best-fit value of the cool temperature in 
the 2T fit. 

From the fit we find that, based on the chi-square (x^) val- 
ues, we cannot significantly differentiate between the wdem and 
the two-temperature model. The wdem model that we fitted had 
the same number of free parameters (upper temperature cutoff, 
lower temperature cutoff, power-law slope, normalization and 
element abundances) as the two-temperature fit (two tempera- 
tures, two normalizations and element abundances) and the dif- 
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Fig. 5. Fraction of the emission measure of the gas below 1.5 
keV to the total emission measure from the wdem fit. 90 cm 
radio contours are overplotted in white. 



ference between the;^^^ values is between -20 and 10 points for 
a typical number of degrees of freedom per bin of 1050-1 100. 

4.3. Discussion of tine temperature structure and 
physical feasibility of the different models 

Although based on the improvement of the values we can- 
not determine which model provides a better description of the 
data, there are several arguments in favor of the continuous 
temperature distribution. 

Firstly, during the rise of the radio lobes, gas at several radii 
should get entrained, mixed, and the uplifted gas should expand 
and cool adiabatically as it moves into regions of lower pres- 
sure. Alternatively, the uplifted gas could be heated up by e.g. 
heat conduction from the surrounding hot phase or by energy 
dissipation from turbulent motions. In any case, it is very un- 
likely that the heating and cooling should exactly balance each 
other to yield the spatially constant temperature of the cool 
component that we observe in the two-temperature fits (both 
MEKAL and APEC). More probably, the constancy of the cool 
temperature could be an artifact of approximating a more com- 
plex temperature structure with a two-temperature model. In 
order to test this, we simulated a series of wdem spectra with 
different slopes, diff'erent high-temperature cutoffs, and a lower 
temperature cutoff' of 0.25 x T^ax, which is a typical value we 
find in our fits, and fitted these simulated spectra with a two- 
temperature MEKAL model. We find indeed, as plotted in Fig. 
[8] that the best fit temperature of the cool component in the 2T 
fit is constant at around 1 keV for a wide range of slopes and 
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Fig. 6. Map of the inverse slope a from the wdem fit. Higher 
values in the map reffect regions where the emission measure 
distribution is ffatter. 
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Fig. 7. Map of the lower temperature cutoff" from the wdem fit. 
Emission from gas at as little as 0.6 keV is detected. Only the 
bins are shown where the fraction of gas below 1.5 keV was 
more than 3cr significant. 
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Fig. 8. Results of fitting various simulated wdem spectra with 
a two-temperature model. The best fit cool temperature in the 
2T model remains constant at around 1 keV for a wide range of 
input parameters in the simulated spectra. The different curves 
are for different values of the upper temperature cutoff, Tj^^x- 



high-temperature cutoffs. Thus the fact that the cool tempera- 
ture in the 2T fit is so constant is very likely an indication that 
the 2T fit is only an approximation of a more complex contin- 
uous temperature distribution. We note that very similar results 
are also found for lower temperature cutoffs of between ~ 0.2- 

0.5 XTniax- 

Secondly, Werner et al. ( 2006a[ l detect Fexvii lines in the 
RGS spectra of M87, lines which are emitted between 0. 14 and 
0.85 keV ( Arnaud & Raymond||1992j l. The presence of these 
lines is inconsistent with the results of the two-temperature fits 
(both MEKAL and APEC) and can only be accounted for with 
a multi -temperature model which predicts emission from gas 
with temperatures of as low as 0.6 keV. Moreover, the RGS 
results show the presence of gas at ~ 0.7 keV out to 3.5' in the 
SW arm, which also speaks in favor of the wdem model. 

Our data thus points towards the conclusion that a contin- 
uous differential emission measure model for the gas in the 
X-ray arms is likely to be a more physical and more realistic 
description than a two-temperature plasma approximation. In 
most of the following sections we will thus primarily focus on 
the results of the wdem fit. 



4.4. Discussion of tine implications of 

multi-temperature structure for thermal conduction 

A gas parcel with volume V and number density n, in thermal 
contact over an area A with an ambient gas with a temperature 
which is different by AT, will reach thermal equihbrium with 



this ambient gas through heat conduction on a time scale that 
can be estimated as (adapted from ,Bohringer & Fabian|1989| l: 



The 



AH 
qA 



InksAT V 
fKoT^'^VT A 



(4) 



where AH is the enthalpy which must be transfered to the 
gas parcel so that it reaches thermal equilibrium with the 



ambient gas, q is the heat flux according to Spitzer (19621, 
kq is the Spitzer heat conduction coefficient, kq ~ 6 X 
lO^^erg cm"' s"' K"^'^, and / is a coefficient allowing for the 
reduction in the heat conduction due to the effects of the mag- 
netic field. For the case of a continuous multi-temperature dis- 
tribution as the one described by the wdem model, the calcu- 
lation of the heat conduction flux and hence the heating time 
scale is generally more complex. We can however make sim- 
plifying assumptions which will lead to upper limits on the al- 
lowed heat conduction such that the cool gas detected survives 
for over ~ 10^ years, which is the estimated age of the inner 



radio-lobes (e.g. Forman et al. 2007 1. The temperature gradient, 
and thereby also the heat flux, is smallest if we assume that the 
gas is stratified in each spatial bin such that gas at the coolest 
temperatures is furthest away from gas at the hottest tempera- 
tures. Moreover, V/A is largest assuming a spherical geometry. 
Combining these two assumptions and the quoted age of the 
radio lobes we arrive at an upper limit for the allowed /. 

We aim to rewrite Eqn [4] which holds for two gas phases 
in thermal contact, for the case of the multi-temperature distri- 
bution described by the wdem model. For this we will calculate 
how much enthalpy AH must be transferred to the gas con- 
tained in the wdem model below a chosen in order to heat 
this gas up to T*. This is given by: 



AH 



2 JcT„,„, 



(«, + ndkniT, - T) dV. 



(5) 



Note that according to the continuous distribution of the wdem 
model, r, can be arbitrarily close to cTj^^^^^, although it is proba- 
bly not physically feasible for the separation to be infinitesimal. 
Given that dY - (dY/dT)dT = UeiiidV and assuming pressure 
equilibrium in the form (n^ + n,) T - const, it follows that 



dV oc T-^^'^dT. 



(6) 



The proportionality constants can be calculated assuming a 
total volume for each bin (see Sect. |6.3| l. The volume V in 
Eqn |4] can then be calculated integrating dV from cT^niax to 
Tf Assuming the gas below to be inside a sphere, one 
can then calculate the area A corresponding to this V. Finally, 
Vr - dT/dr = dT/dV ■ dV/dr can be evaluated at r» using the 
assumption of spherical geometry and Eqn|6] This allows the 
calculation of t/,c, and the requirement that ti,c(T) < 10^ years 
for various r» can help us place upper limits on /. However, 
we cannot account for possible mechanisms that would change 
the temperature cutoffs or the slope of the emission measure 
distribution since the time it was produced. Such mechanisms 
include mixing and the fact that gas at very low or very high 
temperatures, which could have been present at the beginning, 
has already been heated/cooled and is no longer detected. 
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Fig. 9. Factor by which the Spitzer heat conduction must be 
suppressed (this corresponds to 1 // in Eqn.|4]l in order to main- 
tain the cool gas in the multi-temperature structure over at least 
10^ years, as a function of the size corresponding to the volume 
of the cool gas, V(< r»). The wdem inverse slope denoted here 
by p and lower temperature cutoff cr^ax (in keV) of the 5 dif- 
ferent bins used for the calculation are indicated in the legend. 
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Fig. 10. Fe map determined from the wdem model. Beyond the 
expected radial gradient, one clearly sees the enhancement in 
Fe abundance in the arms, especially within the E radio lobe. 
The half-light radius of M87 is marked with a black circle, the 
90 cm radio contours are over-plotted in white. 



We plot in Fig. |9]the conductivity suppression factor 1// 
needed to maintain the cool gas in the multi-temperature struc- 
ture over at least 10^ years against the radius r ~ (3V/4ny^^ 
corresponding to V{< T,). For the calculation we used T, - 
cT^max in the range from 5 x 10^-5 x 10^ K. We calculated the 
values plotted in the figure using results from 5 different bins in 
the E and SW arms spanning values of a between 0.5 and 2. All 
our bins had similar low-temperature cutoffs of cr^ax ~0.64- 
0.76 keV. We find that, for typical blob sizes of 0. 1 kpc, similar 
to the widths of the cool filaments seen in the Chandra im- 
age (the filaments do not appear in images produced in energy 
bands above 2 keV, Forma n et al.||2007) l, a suppression factor 
on the order of ~30-100 is required. 

5. The spatial distribution of metals 

5.1. Abundance maps 

The spectra in each spatial bin provide enough statistics not 
only to test the multi-temperature structure of the gas but also 
to determine the abundances of several heavy elements and to 
analyze their spatial distribution and the influence of different 
multi-temperature models on the abundance determination. In 
the spectral fits we left the abundances of O, Mg, Si, S, Ar, 
Ca, Fe and Ni free. The Ne abundance was fixed to the Mg 
abundance, since the Ne line falls within the Fe-L line complex 
and thus the Ne abundance cannot be accurately determined 
with the spectral resolution of the EPIC detectors. 



Because of the good statistics around the Fe-L line com- 
plex, the Fe abundance is usually determined with the lowest 
statistical errors. We present a map of the Fe abundance de- 



termined from the wdem model in Fig. 10 As expected from 
previous work (e.g. ' Matsushita et al.|2003| , a radial gradient is 
apparent. Beyond this radial trend however, deviations associ- 
ated with the inner radio lobes can be seen as clear evidence of 
the influence of the AGN on the spatial distribution and trans- 
port of metals. The best region to illustrate this is the eastern 
lobe, which from the "mushroom-cloud" shape in the radio im- 
age appears as a more undisturbed system than the SW lobe. 
Here, one clearly sees a stem-shaped enhancement in the Fe 
abundance which coincides very well with the rising plasma 
bubble. The SW lobe also presents a higher metallicity with re- 
spect to the average radial profile, but this enhancement seems 
a bit less concentrated than in the E lobe and is more difficult 
to see from the map in Fig. [TO] A more detailed quantitative 



analysis will be presented in Sect. 6. 1 



Since the Fe-L complex is very sensitive to the temperature 
of the gas, one might question if and how the determined Fe 
abundance systematically depends on the model involved. To 



answer this, we plot in Fig. 1 1 the Fe abundance determined 



from the 2T and the wdem models. Reassuringly, the corre- 
lation between the two sets of fit values is very tight, which 
allows us to be confident in the measured abundance trends. 
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fe (2T fit) 

Fig. 11. Fe abundance determination with two different multi- 
temperature models. The red line represents the function 
Fe{wdem)=Fe(2T). It is easily seen that the abundance deter- 
minations from the two models agree very well, with a very 
low scatter. For clarity, only the points are plotted where the 
emission measure of the cool component in the 2T fit and the 
fraction of cool gas in the wdem fit were more than 3cr signif- 
icant. We note, however, that the correlation is equally strong 
also for all the other points in the maps. 



5.2. Metal abundance patterns of different elements in 
and outside the multi-temperature regions 

The spatial distributions of other elements, especially S and Si 
which have the next best determined abundances, are very sim- 
ilar to that of Fe. Thus, we chose not to include 2D maps for 
these elements but rather to plot their abundance trends relative 
to Fe. We find that, while the overall abundance in the arms is 
higher, the Fe, Si, S and O abundances correlate well both in- 
side and outside the radio lobe regions. The scatter in the O/Fe 
relation is larger than for Si/Fe and the O/Fe slope is much 
shallower, but we do see a small increase in the O abundance 
for higher Fe values. Therefore the cool gas must be enriched 
also with oxygen, for example through stellar mass loss in the 



galaxy. Fig. 12 shows a plot of the Si and O abundances against 
Fe, both for the bins where the cool gas fraction was significant 
(marked with special symbols) and for the rest of the bins in 
the fit. Since the Si and S abundances are very similar, we do 
not add the S/Fe data points to the plot for legibility. Both for 
Si/Fe and O/Fe we see that the trend in the multi-temperature 
regions coincides well with the trend outside the influence of 



gas-uplift by the radio lobes (note that Gastaldello & Molendi 
|2002| also find a constant ratio of O/Fe with radius in the inner 
9' of M87). To quantify this, we fitted a line to the O vs. Fe, Si 
vs. Fe and S vs. Fe data points in and outside the arm regions. 
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Fig. 12. Oxygen and silicon vs. iron abundances in and outside 
the arm regions {wdem fit). Both the Si vs. Fe and O vs. Fe 
in the arms (marked with special symbols, see graph legend) 
follow the same trends, respectively, as the Si vs. Fe and O 
vs. Fe data points for bins outside the arms. Thus, the same 
abundance patterns are seen in all environments. 



The best fit slopes are given in Table [T| These slopes are ex- 
pected to become different at low values of the Fe abundance, 
where the contribution by SNela goes to zero and the O/Fe and 
Si/Fe ratios are determined by the pre-enrichment patterns of 
SNell (for detailed discussion of this see, e.g., Matsushita et al. 



2003 1. With the bins included in our plot in Fig. 12 however, 
we do not seem to reach down to low enough Fe abundances to 
see this effect significantly. 

Table 1. Comparison of best-fit metal abundance ratios in and 
outside the multi-temperature regions (wdem fit). 





inside 


outside 




multi-temperature region 


multi-temperature region 


O/Fe 


0.44 ± 0.07 


0.49 ± 0.03 


Si/Fe 


1.02 ±0.06 


1.10 ±0.02 


S/Fe 


1.05± 0.09 


1.18 ±0.03 



As seen in Table [T[ there are indications that Fe is slightly 
more abundant in the multi-temperature regions (O/Fe, Si/Fe 
and S/Fe are all smaller than in the single-temperature regions), 
which would imply a more important relative contribution from 
SNela, but none of these differences is more than 2cr signifi- 
cant. Because of this and since the O vs. Fe and Si vs. Fe data 
points in the multi-temperature regions overlap well with the 
trend outside the influence of gas-uplift by the radio lobes, we 
conclude that the AGN, at least at this advanced stage of the 
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evolution of the galaxy, enriches the gas halo without altering 
the relative abundance patterns. 

6. The properties of the cool gas 

6.1. The correlation between the amount of uplifted 
cool gas and metallicity 

Several authors have discussed the AGN/ICM interaction as a 
driver for metal-transport from the central galaxy into the ICM 
(for examples of recent work, see e.g. Rebusco et al. 2006[ 
Roediger et al. 2007 1. Comparing Fig. 10 to [4] or [5] a possi- 



ble mechanism by which the AGN influences the distribution 
of the metals in the gas halo is readily apparent: if the cool 
gas from the center of the galaxy is richer in metals, then the 
uplift of this gas by the radio bubbles will also enable the trans- 
port of the metals towards the outskirts of the cluster. Indeed, it 
seems that regions of enhanced metallicity in the lobes corre- 
spond well to regions where the fraction of the cool component 
is large. 



To illustrate this correlation, we plot in Fig. 13 the Fe abun- 
dance as a function of the fraction of the cool component in 
the wdem model for the regions where this fraction was more 
than 3cr significant, for the E and SW arm separately. In the 
plot, the Fe abundance was binned in 2%-intervals of the cool 
gas fraction. We calculated the linear correlation coeflicient for 
each arm without binning, taking into account the error bars on 



both axes ( Akritas & Bershady,1996j , and obtained a value of 
0.66 for the E arm and 0.42 for the SW arm. Combining the 
data points for both arms we obtain a correlation coefficient of 
0.48. This confirms a moderate linear dependence of the mea- 
sured Fe abundance on the fraction of the cool component over- 
all. Furthermore, binning the Fe abundance in 2%-intervals of 
the fraction of the cool component and recalculating the lin- 
ear correlation coefficients, we obtain a value of 0.86 for each 
of the two arms, indicating a very clear correlation. The lin- 
ear fits were performed taking into account the errors on both 
axes. Extrapolating the best-fit line equations to 100% cool- 
component yields an average iron abundance of 1.88 ±0.18 and 
1.93 + 0.12 solar for the S W and E arm respectively (quoted er- 
rors are at the Icr level). 

We adapted this procedure also for the results of the two- 
temperature fits and obtain comparable results for the mean 
iron abundance of the cool gas: 2.13 ± 0.54 for the SW arm 
and 2.38 + 0.29 for the E arm. However, the correlation is not 
as strong as in the case of the wdem model. The coiTelation 
coefficients after binning the data in 2%-wide bins of the cool 
gas fraction are only 0.63 for the E arm and 0.34 for the SW 
arm. The fact that the Fe abundance coiTelates much better with 
the fraction of cool gas determined from the wdem model than 
with the fraction determined from the 2T model may be an ad- 
ditional indication that a continuous temperature distribution is 
a better description of the temperature structure in the regions 
associated with the radio lobes. 

To check to which extent the coiTelation presented in Fig. 
13 may be due to similar radial trends of the cool gas frac- 
tion and Fe abundance, we furthermore fitted the radial Fe and 
/cool profiles from the wdem results using non-parametric, lo- 
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Fig. 13. Averaged Fe abundance in 2%-wide bins of the frac- 
tion of the cool component and linear fit for each arm. The lin- 
ear correlation coefficients are 0.86 for each of the two arms, 
indicating a very clear correlation. 

cally weighted, linear regression smoothing, and divided the 
Fe and /moI maps by the resulting radial models. This should 
remove any trends due to common radial dependences. We plot 



in Fig. 14 the map of Fe deviations from a smooth radial pro- 
file, Fe/ < Fe >, and over-plot contours of the /cool deviations 
from radial symmetry, /cool/ < /cool >■ The figure shows a tight 
correlation between regions of relatively higher iron abundance 
and relatively higher fraction of cool gas with respect to the ra- 
dial average. The only exception to this is a high-abundance 
region to the SE, which does not coiTespond to a high cool-gas 
fraction but where |Simionescu et al.| ( |2007[ ) have already dis- 
covered the presence of a cold front associated with a metallic- 
ity jump. The relatively low-abundance region opposite the SE 
feature is probably due to the fact that the higher abundances 
over a relatively large azimuthal range in the cold front pull 
up the average at that particular radius, so dividing "normal" 
values at that radius by this average will result in a depres- 
sion. We calculated the linear correlation coefficients between 
Fe/ < Fe > and /cool/ < /cool > and obtain 0.70 for the E 
and 0.72 for the SW arms (unbinned). The fact that these co- 
efficients are even stronger than for the simple Fe vs /cod cor- 
relation indicates that the underlying radial dependence of the 
abundance of the hot halo may be responsible for a large part 



of the scatter in Fig. 13 We will investigate this further using 
spectral simulations in the next section. 

6.2. The metallicity of the cool gas 

To explain the correlation between the amount of uplifted cool 
gas and metallicity presented in the previous section we can 
begin by assuming that if we fit a multiphase gas for which the 
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Fig. 14. Map of Fe deviations from a smooth radial profile, 
Fe/ < Fe >. Contours of the cool gas fraction deviations from 
radial symmetry are over-plotted. It is clear that the correlation 
between Fe and /cooi is in most part not due to common ra- 
dial dependences. The high-abundance feature to the SB corre- 



sponds to the cold front presented in Simionescu et al. (2007 1. 



Fe abundance is some function of the temperature with a spec- 
tral model with a single metallicity, the resulting Fe abundance 
from the spectral fit can be (neglecting for now e.g. the depen- 
dence of the Fe line strength on temperature) approximated by 
an emission-weighted average: 



X^"" Zpe(T)Y(T)dT 

^ mux 

rT^ Y(T)dT 

"-^C i max 

(^Fctuoi ) /cool + (^Feh„, ) /hot 
(•^Fcho,) + ((ZFecool) ~ (•^Fei,„,))/cool 



(ZFec„ol ) /cool + (Zpeto, > /hot (7) 



(8) 
(9) 



where the quantities denoted by the prime are the results of the 
fitting while the other quantities are the "true" values charac- 
terizing the gas, and it holds that f^^^ + f^^^ - 1 since /J^^j and 
/jj^j denote the fitted relative emission measure contributions 
from the cool and hot gas respectively. The cool gas is, as be- 
fore, defined as gas with a temperature below 1.5 keV, which 
is the lowest value found in the deprojected temperature profile 
outside the X-ray arms (Matsushita et al."2002 ). 

This means that the measured Fe abundance Zp^ depends in 
first-approximation linearly on the fraction of the cool, metal- 
rich gas /(,'oop or equivalently on Y^ooi/Y, which is in agreement 
with what we observe. If this approximation holds, then the Fe 
abundances obtained by extrapolation to 100% cool component 
contribution, as calculated above, are the real values of the Fe 
abundance of the cool component. 

To investigate the validity of this simple approximation, 
we simulated spectra combining a wdem component between 
c Tmax and 1 .5 keV with Fe = 2 solar and a second wdem com- 
ponent between 1.5 keV and r,nax with Fe - [0.6,0.8,1.0] so- 
lar. The two wdem components had the same slopes and their 
relative normalizations were chosen such that dY/dT at 1.5 
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Fig. 15. Results of fitting a simulated combination of two wdem 
spectra with a high Fe abundance of the cool gas with a wdem 
model with coupled abundances. Blue, black and red crosses 
plot the fitted (coupled) Fe abundance against the fitted fraction 
of the cool component for various input wdem slopes, various 
upper temperature cutoffs, an input lower temperature cutoff 
of 0.25 r„ax ™d a hot component metallicity of 1.0, 0.8 and 
0.6 solar, respectively. Black circles, black triangles and black 
asterisks plot the same for a hot component metallicity of 0.8 
solar and lower temperature cutoffs of 0.2, 0.3 and 0.4, respec- 
tively. Green squares plot the fitted (coupled) Fe abundance 
against the original (input) cool gas fraction for an input lower 
temperature cutoff of 0.25 T^ax and a hot component metallic- 
ity of 0.8 solar For an ideal fitting of the cool gas fraction, the 
black crosses should overlap with the green squares. 

keV was continuous. We stepped the upper temperature cutoff 
Tmax between [2.2,2.6,3.0] keV, the lower temperature cutoff 
between [0. 2,0.25, 0.3,0.4]rniax and the slope between 0.2 and 
1.2 in intervals of 0.2. This parameter space is equivalent to 
a range of cool gas contributions from to 60% of the total 
emission measure. We assumed an O/Fe ratio of 0.4 and Si/Fe 
and S/Fe ratios of 1.0 solar (see Table [T} for both of the sim- 
ulated model components. The simulated spectrum was then 
fitted with a single wdem model with a single abundance. 



The results are presented in Fig. 15 Most important to note 
is that from the simulation we do recover, up until quite high 
percentage contributions of the cool component, a clearly linear 
dependence of the determined coupled abundance Zp^ on the 
cool gas fraction /Joop plotted on the X axis. The largest scatter 
is due to the different abundances of the hot gas component, 
while there is practically no scatter due to c and T^ax- However, 
the slope of Zp^ vs. /'^^j is slightly different than the slope of 
Zpg vs. the true /ooi (by a factor of approximately 1.2). 

This discrepancy can be understood in the following way: 
if the cold gas is more metal-rich than the hot gas and if the 
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spectral model forces the abundances of all gas phases to be 
the same, then the cool gas in the fitted model will have a too 
low abundance compared to the real value (and in turn the hot 
gas will have a too high abundance). To obtain the best fitting 
spectral shape of the Fe-L complex thus, which is equivalent 
to obtaining the correct contribution of emission in the FeL en- 
ergy range from the cold gas, the fit will have to raise the nor- 
malization of the cool component to compensate for the fact 
that the fitted abundance of the cool component is forced to be 
lower than in reality. This is why coupling the abundances of 
the different phases leads to an overestimation of the normal- 
ization of the cool component, which is reflected in our simu- 
lations by a decrease of the Zp^ vs. /^'^^j (fitted) slope compared 
to the vs. /cool (real) slope. This affects our estimates of the 
normalization of the cool gas, as a consequence also our esti- 
mates of the mass of the cool gas, and leads to underestimating 
the metallicity of the cool component if we simply extrapolate 
the best-fit Z' vs f , line equations to f' , = 1 . If we correct 

re ■'cool ^ ■'cool 

the slopes found in the previous section for the wdem model by 
the factor described above we find Fe abundances for the cool 
component as summarized in Table |2] 

We note that similar simulations as described above were 
performed also for the two-temperature model, recovering 
again the expected linear dependence of Zp^ on the cool gas 
fraction /^'oop but with a larger scatter due to the assumed metal- 
licities of the hot component. Given the much weaker corre- 
lation between Zp^ and /J^^j in the 2T model and the higher 
uncertainty in the determination of the observed Zp^ vs. /J^^j 
slope, as well as considering our arguments against the two- 
temperature description of the plasma the X-ray arms, we how- 
ever choose not to further discuss the 2T results in quantita- 
tively determining (Zpe^^^^,). The assumed thermal model may 
have an important influence on the determination of the value 
for (ZFe„„|). However, more complicated models cannot be fit- 
ted to the available data, since any additional free parameters 
would be unconstrained. In view of the current lack of a clear 
theoretical prediction regarding the exact thermal structure of 
the uplifted gas, the wdem model is at the moment the sim- 
plest and most flexible model which serves as a good first- 
approximation to a wide range of possible emission measure 
distribution functions. Hence, we will in the following rely on 
the wdem results for interpretation. 



Table 2. Estimates of the average iron abundance of the cool 
component relative to the proto-solar values of Lodders ( 2003 | l. 





E arm 


SW arm 


2T 


2.38 ± 0.29 


2.13 + 0.54 


wdem uncorrected 


1.93 ± 0.12 


1.88 ±0.18 


wdem corrected 


2.16 ±0.14 


2.11 ±0.22 



Considering the results shown in Table |2] we will for the 
rest of this work use an average iron abundance of 2.2 solar. 
This impUes, if we use the abundance ratios presented in Table 
[l] Si and S abundances also around 2.2 solar, and an O abun- 
dance of roughly 0.9 solar for the cool gas. Interestingly, the 



determined Fe abundance for the cool gas is in good agreement 
with the value of 2.3 solar obtained by [Matsushita et aL ( 2007| l 
in the Centaurus cluster, where presumably an undisturbed ac- 
cumulation of metals over a longer time-scale enabled the cen- 
tral metallicity to reach higher values than in the Virgo cluster 
center. The difference in the central abundances of M87 and 
Centaurus can be thus explained by AGN-induced transport of 
heavy elements out of the central region of M87. 

6.3. The mass of the cool gas 



The emission measure is defined as 
Y - I ni.n\dV ~ ngtiiV ^ feinlv 



(10) 



where «e and rii are the electron and ion densities, /e; is the 
electron to ion number ratio (on the order of 1 .2), and V is the 
volume of the emitting region. From this, one can estimate the 
mass of gas according to (where (m,) is the average ion mass 
per electron in the cluster gas and nip is the proton mass): 



(11) 



This requires assumptions about the projected geometry. To 
calculate the mass of cool gas in the arms, we compare two 
alternatives for determining the depth of each bin along the line 
of sight (LOS). 

6.3.1 . Assumptions about the projected geometry of 
the cool component 

The first and simplest method is to assume a constant LOS 
depth for the cool gas component equal to the average width 
of the radio arms in the plane of the sky. For this, we choose 
a LOS depth of 200" (15.5 kpc). The second possibility is to 
define the center of the gas halo, determine the minimum and 
maximum radii of each bin with respect to this center and as- 
sume that the longest LOS distance is equal to the longest cir- 
cle chord that can be drawn between the circle of minimum 
radius and the circle of maximum radius, in other words that 
only the gas between the sphere with the minimum radius and 
the sphere with the maximum radius contribute to the emission 
( [Henry et al.|2004[|Mahdavi et al.|2005| ). This yields, where S 
is the area of the region in the plane of the sky. 



L = 2 



.2 

max 



mm' 



and 



V^2SL/3. 



(12) 



We used both of these geometries combined with the inte- 
grated emission measure y^ooi of the gas below 1.5 keV from 
the wdem fit to calculate different estimates of the mass of the 
uplifted cool gas. In the calculation we included only the bins 
where the fraction of gas below 1.5 keV was more than 3cr 
significant. Both methods show similar spatial distributions of 
the mass of cool gas. The quantitative results for the total mass 
of the cool gas are summarized in Table |3] under the column 
heading "direct". For comparison, also the masses of the cool 
gas determined from the two-temperature model with the two 
assumed geometries are shown and are in good agreement with 
the wdem results. 
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6.3.2. Filling factor 



For the two volume estimates considered above we additionally 
must take into account the fact that the cool component may not 
fill the entire volume of the bin, but may be concentrated in thin 
filaments ( [Forman et al.|[2007 i which occupy only a fraction 
(so-called filling factor) of the assumed volume corresponding 
to each bin. It is very difficult to determine this filling factor 
from the data, especially since the PSF of XMM-Newton is too 
large to resolve these filaments. 

The easiest method to estimate the filling factor is to as- 
sume pressure equilibrium between the different temperature 
phases present in the spectral models, which for the wdem 
model implies SV oc f^+Ua (ggct. 4.4 1, yielding 



(13) 



The results for the new mass estimates using the filling fac- 
tor corrections according to these calculations are presented in 
Table[3]under the column heading "ff". We used only the vol- 
ume estimate given in Eqn [T2j since we cannot assume that 
the hot component also has a constant line of sight distance of 
only 200". These values are however only lower limits to the 
true mass of the cool gas because the filaments are most likely 
supported by extra pressure from the magnetic field. If the mag- 
netic field is higher inside the cool component than in the sur- 
rounding hot component, then the magnetic pressure adds to 
the thermal pressure of the cool gas. This reduces the pressure 
required to maintain hydrostatic equilibrium with the hot sur- 
rounding gas, which allows a lower density of the cool gas and 
a larger filling factor. Note that we have ignored the effects of 
magnetic tension which can contain the cool gas and work the 
other way. 

To get a more physical estimate of the filling factors we 
therefore also reevaluated them including the contribution to 
the pressure of the cold phase of a uniform magnetic field with 
a strength of 10 fiG, which is a typical average value in the inner 
radio lobes obtained by jOwen et al. ( 2000 
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Fig. 16. Map of the mass of gas below 1 .5 keV in the wdem 
model assuming a geometry with a line-of-sight distance for 



each bin given by Eqn 12 The map has been corrected for the 
filling factor of the cool gas assuming pressure support for the 



cool gas from a lO/uG magnetic field (Sect. 6.3.2 1 and for the 



effects of coupling the abundances of the different phases (Sect. 



6.3.3 1. The units are solar masses per square kiloparsec. Radio 



contours at 90 cm are overplotted in white. 



2.2 solar and find that this factor changes by only 5%). We 
therefore further correct the mass of the cool gas decreasing it 
by this factor, as shown in Table |2] under the column heading 
"ffBlO Zcorr". A map of the mass of cool gas including all the 



using the radio data, corrections described above is shown in Fig. 16 



The results are likewise shown in Table 3] under the column 
heading "ff BIO". We note that the lower limits on the magnetic 
field strength of 1 yuG obtained in this work (Sect. [7]) would 
result in a magnetic pressure several orders of magnitude below 
the thermal pressure, therefore we did not include this in our 



Table 3. Total mass of the cool gas (in 10^ solar masses), in- 
cluding all the different corrections discussed in Sect. 6.3 



calculations but preferred the higher values obtained by jOwen 
|etaL] ( [2000l i. 

6.3.3. Corrections due to coupling the abundances of 
different temperature phases in the spectral 
models 



LOS 


model 


direct 


ff 


ff 
BIO 


ffBlO 
Zcorr 


const 


2T 


2.3 








Eqn 


12 




2T 


2.4 


0.39 


0.41 




const 


wdem 


2.2 








Eqn 


12 




wdem 


2.1 


0.54 


0.60 


0.50 



We proved in the Sect. 6.2 that the normalization of the cool 
gas is overestimated by using a model with a single abundance 
to describe a combination of a metal-rich component and rela- 
tively poorer hot ambient plasma. In the case of a wdem model, 
this overestimation amounts to a factor of 1.2 if the average 
metallicity of the cool component is 2 solar (we redid the sim- 
ulations for an average metallicity of the cool component of 



In the following, we will use the estimate of 5 x IO^Mq 
for the mass of the cool gas, which includes all the correc- 
tions discussed in this section and is determined from the wdem 
model which we believe to describe the data better than the 
two-temperature model. We note that the cooling time of the 
cool gas in the multi-temperature regions, computed assuming 
isobaric cooling and using an average value of the temperature 
for the multiphase distribution, is longer than the ages of the 
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radio lobes in all bins (between 2 x 10** and 2 x 10^ years). 
Therefore our mass estimate is not affected by gas which has 
cooled out of the X-ray domain since the uplift. 

6.4. Discussion of the origin and uplift of the cool gas 

We have presented so far in this work many arguments in favor 



of the scenario proposed by Churazov et al. (2001 1, in which 



the cold gas in the X-ray arms is a result of gas uplift from the 
central regions of the galaxy through buoyant radio bubbles. 
In this section we will try to bring the observational details to- 
gether towards a consistent picture of the energetics and chem- 
ical composition of the M87 X-ray arms and explain how this 
relates to the radio activity. 

6.4.1. Energetics 

We first use the mass estimates derived in the previous section 
to compute the gravitational energy needed to uplift the gas: 



U = 



GM,,,iM(R) 
R 



(14) 



For M(R) we use the integrated total mass profile of Matsushita 



et al. ( 2002 1. We use the approximation that all the gas with 
total mass Mcooi is uplifted from R=0 and U=0 to a mean radius 



i:M(i)R(i) 



(15) 



where i are all the bins where the cool component was more 
than 3cr significant. We calculated Rmean in this way for all the 
different models and volume estimates presented in the pre- 
vious section and obtain consistent values between 15 and 20 
kpc. Assuming a typical mass of the uplifted gas of 0.5x lO^M© 
(see Table|3]l and ^mean = 20 kpc, we obtain a required gravita- 
tional energy of 4.33 x 10^^ ergs. If the AGN outburst that cre- 
ated the inner radio lobes occurred 10^ years ago (e.g. Forman] 
et al.|2007 1, this gravitational energy corresponds to a required 
power of 1.37x10"*^ erg/s, which is a sizable fraction of the total 
jet power estimated by Owen et al. (2000) at a few xlO'*'* ergs/s. 
The power needed to upli ft the cool gas is m uch larger than the 
total radio power cited by 'Owen et al. (2000) at only 9.6 x 10"*' 



erg/s, consistent with the picture that the bulk of the jet energy 
is deposited in the gas (in this case, a large part thereof in the 
form of potential energy), while only a small fraction is leaked 
to radio emission. 

Thus, from energetic considerations, the model of the cool 
gas being uplifted by the buoyant radio lobes is entirely con- 
sistent with the X-ray and radio observations. This provides a 
clear mechanism for the transport of the cool gas out of the cen- 
ter of the gas halo, together with the metals it has been enriched 
with. 

6.4.2. Possible sources of cool gas 

We have chosen to define the cool gas in the wdem model as 
gas having a temperature below 1.5 keV, the smallest value 
in the deprojected profile of the hot M87 halo calculated by 



Matsushita et al. ( 2002| l. There are three possible ways to pro- 
duce gas cooler than 1.5 keV. The most straight-forward option 
is to adiabatically cool gas from the center of the M87 halo by 
uplifting it into regions of lower pressure at larger radii. Using 
the pressure profile of [Simionescu et aL ( 2007[ ) and assuming 
an adiabatic index of 5/3, gas at 1.5 keV uplifted from a ra- 
dius of 1 to 10 kpc would cool from 1.5 to 0.6 keV and so 
could reproduce, in principle, the temperatures that we observe. 
However, if adiabatic cooling is the only mechanism to reduce 
the temperature of the gas, we would (neglecting heating pro- 
cesses) expect to see a decreasing trend of crmax with radius, 
which is not the case. More importantly, the average Fe abun- 
dance of 2.2 solar for the cool gas, determined in Sect. 6.2 is 
higher than the value observed in the center of M87 (~ 1 .6 solar. 



Matsushita et aL|2003| l. Thus, the cool gas could not originate 
from adiabatically cooled halo gas alone, but must contain a 
second, metal-rich component. 

Another important source of cool gas are the stellar winds. 



Assuming the stellar mass loss rate proposed by Ciotti et al. 
( fT99T] l, 



MAt) X -1.5 X 10-^^ Lb t-^^ . 



(16) 



where Lb is the blue band luminosity in units of L^o and fi5 is 
the age of the stellar population in units of 15 Gyrs. We take 
the age of the stellar population to be 10 Gyrs and Lb to be 
1O"Lb0, corresponding to the absolute B magnitude of M87 
of -22.14 cited by Peletier et al.| ( |1990) . This gives us a stel- 
lar mass loss rate of 2.54 Mg/yr. If we further assume that the 
material lost from stars thermalizes according to the M87 stel- 
lar velocity dispersion of 350 km/s ( Angione et al.|1980 i, the 
temperature of this gas would be 0.64 keV. Interestingly, this is 
also very close to the lower temperature cutoffs found using the 
wdem spectral model. Beside being a source of cool gas, stellar 
winds also contribute to the metal-enrichment, as will be dis- 
cussed in the next section, and are an important candidate for 
the second, metal-rich component of the uplifted gas. 

The third and last possible source of cool gas is the cooling 
flow. jMatsushita et al.| ( |2002) l find that within 4' a total mass of 
13.8 Mo/yr can cool down to as low as 1.4 keV and 0.63 Mo/yr 
cool down to 0.1 keV. Based on the same metallicity argument 
as above, the cool gas could not originate from cooling-flow 
gas alone since this would lead to its metallicity being lower 
than measured. Using the gas metallicity and temperature, it is 
not possible to distinguish adiabatically and radiatively cooled 
ICM as different components of the uplifted gas. 



6.4.3. Metal enrichment of the cool gas 

Metals can be either ejected from the stars by the stellar winds 
or produced by supemovae. We neglect the metal enrichment 
by type II supemovae, since M87 has an old population where 
the expected SNII rate is very low. [Rafferty et aL ( 2006[ l, for 
example, find an upper limit on the star formation rate (SFR) 
for M87 of 0.081 Mo/yr using far-infrared data. Assuming a 
Salpeter initial mass function (IMF, [Sa"lpeter| 1 9531 ) for the stars 
forming between 0.01 and 50 solar masses, the number of stars 
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formed from this cooling gas that would be more massive than 
10 solar masses and could become SNll progenitors is 



abundance in the multi-temperature regions vs. the fraction of 



A^(> lOMo) = 0.08 



r50 . , , 



-2.35 



dM 



M M-^-^^dM 



0.017/century (17) 



This upper limit on the rate of SNell is roughly two orders of 
magnitude below the expected rate of SNela in elliptical galax- 
ies ( C^peniro]etal]T999 see calculation below). 

This leaves stellar winds and type la supernovae (SN la) 
as the main sources of enrichment of the cool gas. Following 



Bohringer et al. ( 2004 1; Rebusco et al. ( 2006 1, the rates of iron 
injection by stellar mass loss is: 



^Fe,. - TFe X M,(f), with 



7 X X 

•^Fe.t ^ ~ T >^ 



A^F 



(18) 
(19) 



where M»(f) is the stellar mass loss rate calculated in Eqn. 16 
7Fe is the mean iron mass fraction in the stellar winds of an 
evolved stellar population, Zpe,* is the assumed abundance of 
the stellar winds in solar units, (mpe) ~ 56 is the mean weight 
of an iron isotope with respect to hydrogen, (m,o) is the mean 
weight of an ion with respect to hydrogen in the Sun (in the 
solar photosphere), and (A^Fe/ 2/M)Lodd ~ 3.16 x 10"'' is the 
relative number abundance of Fe ions to the sum over all ions in 
the plasma, as determined from Lodders (2003). We assume, as 
above, a galactic age of 10 Gyr Kobayashi & Arimoto_( 1999) 1 
derive stellar metallicities in M87 ranging from *2 solar in the 
center to ^1 solar at the galaxy's half-light radius (2') using the 
Mg2 index. We will therefore take Zpe,* ~ 1-6. The rate of iron 
production by SN la can be written as: 



M, 



Ve,Ia = T/Fe X/?/fl,with 

R,a = 0.18 ± 0.06 (100yr)-'(lO'%Q) 



(20) 
(21) 



where Jjfe = Q.19Mq is the iron yield per SNla assuming the 
WDD2 model of ( [Iwamoto et al.|[T999) l and Rja is the present 
supern ova rate in elliptical galaxies according to 'Cappellaro 
et al.| ( [T999) . For M87, Lb ~ 10"Lbo, as calculated above, 
which yields Rja ~ 1.8 + 0.6 SNela per century. We assume 
that the iron production by both SN la and stellar winds is con- 
stant during the time of enrichment, which as we will show 
below is much shorter that the Hubble time. 

We wish to enrich ICM gas to an average Fe abundance of 
2.2 solar. For this, we must assume an initial metallicity Z, of 
the ICM at the center of M87 immediately following the next 
to last AGN outburst which modified the spatial distribution of 
metal abundances. After the next to last AGN outburst, stel- 
lar mass loss and supernovae enriched this central gas until it 
reached 2.2 solar and was uplifted by the last AGN outburst, 
the effect of which we observe today. Ideally we should take 
the central metallicity immediately after the last outburst and 
assume it is representative also for the metallicity after the next 
to last outburst. We however do not know this value. A sensi- 
ble lower estimate is to assume Z, as the average metallicity of 
the entire hot halo within the region of the arms. This can be 
calculated as the best-fit y-intercept for the linear fit of the Fe 



cool gas (Fig. 13 1 and is approximately 0.7 solar. Alternatively, 
as an upper limit to Zi, we can assume the current value at the 



center of M87 determined by Matsushita et al. ( 2003 1 to be 1 .6 
solar. This is likely to be an upper limit because, since the time 
of the last outburst, the central region has been enriched with 
metals beyond the metallicity that it had immediately after the 
outburst. With these two values for Zi, we can bracket the en- 
richment time needed to reproduce the observed Fe abundance 
of the cool gas. 

Using a mass production rate of stellar winds of 2.5 Mo/yr 
(Eqn. [T6| and neglecting the mass of supernova remnants (as- 
suming an average supernova mass of 10 Mq and rate Rsn la of 
1.8 per century, we obtain only 0.18 Mg/yr which is much less 
than for the stellar winds), we can write for the Fe enrichment 
of the uplifted material: 



2.2 X 5 ■ lOVo X /Fe,Q = Zi X Mi X /pe.o 

+ (Mpe.Ia + Mpe,.) X T 

Mi + 2.5Mo/yr x r = 5 ■ IO'^Mq 



(22) 



(23) 



where /peo = (mp^/mie) X (A^Pe/ 2/ A^/)Lodd is the mass frac- 
tion of Fe in the Sun, r is the time over which the enrichment 
takes place. Mi is the initial mass of cluster gas to which stellar 
mass loss is added over time t. We solve the system of equa- 



tions above for t and Mi making use also of Eqns. 18 and 20 
Assuming Zi=1.6 solar we obtain t=30 Myr, while for Zi-0.1 
solar we need 60 Myr to enrich the ICM before uplift. For 
Zi=1.6 solar, the mass originating from stellar mass loss rep- 
resents 15% of the total uplifted mass, while for Zi=0.7 solar it 
amounts to 30%. 

We can now try to reproduce the abundances of other el- 
ements such as O, Si and S in the cool gas. Using the Fe/O, 
Fe/Si, Fe/S ratios in Table[T] the O, Si and S abundances of the 
cool gas should be approximately 0.9, 2.2 and 2.2, respectively. 
For an initial Fe abundance of Zi - 0.7 solar and assuming the 
same ratios, we obtain 0.3, 0.7 and 0.7 solar for O, Si and S, 



respectively. The central values of Matsushita et al. ( 2003 j l are 
consistent with 0.6,1.6, and 1.6, respectively. 

Assuming solar abundance ratios for the stellar winds, we 
can use the equivalent of Eqn. 18 and 19 to determine also 
Mo,t. Using Zo,i=0.6 and ignoring contributions by SN la to 
the oxygen budget, the time needed to enrich the cluster gas 
to 0.9 solar oxygen abundance is 60 Myr. Correspondingly for 
Zo,i=0.3, we obtain 92 Myr. A similar procedure can be ap- 
plied for Si and S, including contributions by SN la according 
to the WDD2 model of 0.21 Mq and 0.12 Mg per supernova, 
respectively. For Si we obtain, respectively, 63 and 108 Myr for 
•Zsi,i=1.6 and 0.7, while for S we obtain 56 and 98 Myr for the 
same Zs,i. 

We also computed the required enrichment times using the 
WDDl model of Iwamoto et al. ( \999) which has a higher 
Si/Fe ratio and was found by Finoguenov et al. (2002 1 to pro- 



vide a better fit to the abundance patterns in the center of the 
M87 gas halo. These results, together with the results from the 
WDD2 model, are summarized in Table [4] Oxygen enrichment 
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times for the two models are the same since we neglect con- 
tributions by SN la to the oxygen budget. The WDDl model 
requires more similar enrichment times to produce the differ- 



ent metals, in agreement with the results of Finoguenov et al 
( |2002| l. The obtained enrichment times differ from eachother 
by up to a factor of « 2.1 for the WDD2 and « 1.7 for the 
WDDl models. A possible source of this discrepancy is that 
our assumption of solar abundance ratios for the stellar winds 
is not accurate. If, for the older population of M87, the Si/Fe 
composition of stellar winds would be higher than solar, we 
would obtain more consistent enrichment times to produce all 
three metals. Note also that the enrichment times are in slightly 
better agreement for Z\ - 0.7 solar than for Z\- 1.6 solar, sug- 
gesting that the former choice may be the more appropriate and 
reliable one. 



Table 4. Enrichment times needed to produce the metals seen 
in the cool gas (in Myr) and, in brackets, the fraction of the 
stellar mass loss contribution to the total mass of uplifted gas. 





Zpc/Si/S.i - 


1.6 


Zpc/si/s.i - 0.7 




Zo.i = 0.6 




Zo.i = 0.3 






WDD2 


WDDl 


WDD2 


WDDl 


Fe 


30 (0.15) 


35 (0.18) 


60 (0.30) 


69 (0.35) 


Si 


63 (0.32) 


49 (0.25) 


108 (0.54) 


90 (0.45) 


S 


56 (0.28) 


42 (0.21) 


98 (0.49) 


79 (0.40) 





60 (0.30) 


60 (0.30) 


92 (0.46) 


92 (0.46) 




Fig. 17. Regions for spectral extraction used to test the presence 
of a power-law component (in black) overplotted on the map of 



the gas pressure deviations from radial symmetry (Simionescu 
|et al.|20d7] l. 90cm radio contours are also overplotted in white. 



At the end of Sect. |5.2| we show that the relative abun- 
dances of O/Si/ S/Fe are very similar in and outside the multi- 
temperature regions. This indicates firstly that the metals out- 
side the multi- temperature region, which were (according to 
our model) uplifted from the galaxy by previous AGN out- 
bursts, must have been transported after the last major epoch 
of star formation, when contributions of SNIl to the chemi- 
cal budget of the galactic gas became negligible. Secondly, the 
SN la enrichment must have remained approximately constant 
in time relative to the enrichment by stellar winds. 

We note that |Owen et al.| (|2000 ) estimate the age of the 
outer radio halo between 96 and 150 Myr/P44, where f 44 is the 
total jet power in lO"*"* erg/s. If the total jet power is around 
1-3 xlO''^ erg/s, we would obtain consistent values with the 
enrichment times of ~30-110 Myr obtained above. Thus, it is 
possible that both the outer and the inner radio halos gener- 
ated similar outbursts ~ 3-11x10^ and ~ 10^ years ago respec- 
tively, transporting metal-rich gas outwards, while between the 
two outbursts the metals could accumulate in the center of the 
gas halo. Less violent outbursts which do not have such large 
impact in the distribution of metals in the gas halo could have 
occurred between these two events, as is suggested by the mul- 

( |2007] l. 



7. Constraints on possible non-tliermal emission 



tiple shocks found in the Chandra data by Forman et al 



Note that if some of the uplifted metal-rich gas fell back 
to where it originated between outbursts, the time intei-val be- 
tween two such outbursts could be shorter since one would not 
have to wait for all the extra metals observed in the cool gas to 
be produced from stellar winds and supernovae. 



Simionescu et al. ( 2007] l suggest the presence of non-thermal 
pressure support by relativistic electrons in the X-ray arms of 
M87. We performed a more detailed analysis to check for spec- 
tral signatures of inverse Compton (IC) scattering of the mi- 
crowave background photons on the relativistic electrons in the 
M87 gas halo. For this, we chose 7 annuli between 1 and 9 ar- 
cminutes and divided each of these into 4 sectors containing 
the E aim, the SW arm, and the NW and SE regions between 
the arms. The sectors containing the arms were fitted with 
a two-temperature VMEKAL model and a power-law while 
the NW and SE regions were fitted with a single-temperature 
VMEKAL model plus a power-law. Our extraction regions can 



be seen in Fig. 17 



We find that the fit cannot constrain the power-law index 
in any of the extraction regions, therefore we fixed it to 1.7, 
which is consistent with the number distribution of relativis- 
tic electrons in our Galaxy. With a fixed power-law index, we 
find that the normalization of the power-law component is more 
than 3cr significant everywhere except in the 7th ring of the S W 
arm sector. The flux of the power-law component in all regions 
is a factor of more than 10 above the cosmic X-ray background 
(CXB) found by |De Luca & Molendi] ( |2004) l in the 2-10 keV 
band and a factor of more than 5 above the same CXB model in 
the 7-10 keV band. However, we cannot rule out the possibil- 
ity that the power-law component is significant simply because 
it compensates for calibration or spectral model inaccuracies, 
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therefore we cannot confirm it as a detection of IC emission 
from the M87 halo, but interpret our results rather as an upper 
limit of the IC flux that is consistent with the data. 

For the upper limit on the non-thermal pressure, we use 
( |Lieu et al.|1999l l 



8 X 10^'L42 



2-^i 



p = 


E 

— \ 
3V 


max 


2-fi 

-J 

' min 


max 


-J 

' mill 



where 



and 



2n 



y = y (1 - cos eXR;;,,, - R;J. 



(24) 
(25) 
(26) 



L42 is the luminosity of the non-thermal emission between ymax 
and ymin in units of 10"*^ ergs, is the index of the electron 
number distribution which is related to the photon index of our 
power-law model as ju = 2r - 1 = 2.4, y^in and jmux are the 
minimum and maximum Lorentz factors of the electrons con- 
tributing to the IC emission, V is the volume of a spherical 
shell with half-opening angle 6 and minimum and maximum 
radii 7?n,in and Rmax- We extrapolated our power-law luminosity 
to include IC scattering by electrons with Lorentz factors be- 
tween 7niin - 300 and y^ax - 10000 and find 90% confidence 
upper limits on the non-thermal pressure as presented in Table 
|5] Note that the pressure calculated in this way is only due to 
the relativistic electrons, and that relativistic ions may give rise 
to additional non-thermal pressure up to 30 times higher than 
that of the electrons. We also computed the thermal pressure as 
p = (1 -I- 1 lfi.i)nekT with He given by Eqn[TO]in order to present 
also the relative contribution of the non-thermal pressure to the 
total (thermal H-non- thermal) pressure in Table |5] 

Table 5. 90% confidence level upper limits on the non-thermal 
pressure (in 10"" dyn/cm^ and in brackets as percentage of the 
total pressure). 



^ mean 


NW 


SE 


E 


SW 


(arcmin) 






arm 


arm 


1.38 


7.4 (31.5) 


27.6 (54.4) 


18.2 (40.9) 


19.2 (40.5) 


2.09 


3.2(19.8) 


7.9 (30.1) 


6.7 (25.7) 


8.5 (32.2) 


2.75 


1.7 (14.1) 


2.8(16.5) 


4.5 (22.6) 


2.3 (13.4) 


3.63 


0.85 (9.4) 


1.2(10.3) 


2.1(17.3) 


2.1(15.9) 


4.92 


0.68 (9.9) 


0.50 (5.5) 


1.1 (12.6) 


1.3 (13.6) 


6.55 


0.22 (4.4) 


0.42 (6.1) 


0.53 (7.5) 


0.84(11.3) 


8.25 


0.14 (3.5) 


0.22 (4.1) 


0.31 (5.7) 


0.26 (4.5) 



The upper limits presented in Table |5] are consistent with 
the contribution of the non-thermal pressure of 5-10% sug- 
gested by'Simionescu et al. (2007). However, from these upper 
limits no difference between the arm and ofF-arm regions can 
be seen, suggesting that perhaps at least partly the normaliza- 
tion of the power-law is determined by calibration uncertain- 
ties. The high contribution of the power-law in the inner regions 
is probably due to the central AGN, which was very bright dur- 
ing this observation and could influence the spectrum out to 
larger radii due to the tail of the XMM PSF. 



We also used our upper limits on the power-law flux, as- 
suming it is due to IC emission, to determine the lower limits 
on the magnetic field stength, which can be obtained from the 
equation ( |Sarazin|1986 1 



fr 



2.47 X io-i''r>o«) /4960r, 

Baipi) ^ \ B 



(27) 



where a(ju) and b(/i) are unitless functions of fi only, and for 
our case fi - 2.4, a - 0.086 and b - 7.0, = il-i - l)/2 is the 
X-ray spectral index, Tr is the temperature of the microwave 
background radiation (2.73 K), B is the magnetic field strength, 
ft is the X-ray flux density at frequency Vx in ergs/s/cm^/Hz 
and /, is the radio flux density at frequency v^. We used the 



radio map of |Owen et ar.| ( ,2000) at 90 cm (vr=327 MHz) and 
chose to calculate the X-ray flux at 6 keV by dividing the 5 - 
7 keV power-law flux by the corresponding bandwidth in Hz. 
Our lower limits on the magnetic field strength B are summa- 
rized in Table |6] We could only determine the values for B in 
the first 6 rings (5 for the E arm) because the other regions were 
outside the radio map. 

Table 6. 90% confidence level lower limits on the magnetic 
field strength ifJ.G). The sixth ring in the E arm was outside the 
radio map 



(arcmin) 


NW 


SE 


E 
arm 


SW 
arm 


1.375 


0.53 


0.31 


0.51 


0.65 


2.085 


0.62 


0.48 


0.77 


0.94 


2.750 


0.56 


0.64 


0.96 


1.40 


3.625 


0.60 


0.78 


1.26 


1.09 


4.920 


0.47 


0.82 


0.52 


1.16 


6.545 


0.29 


0.52 




1.20 



Interestingly, we do find more stringent lower limits on the 
magnetic field in the X-ray arms, these being the only regions 
where our data require magnetic fields stronger than 1 /iG. 
However, as these values are only lower limits, we cannot draw 
definitive conclusions based on this fact. 

Assuming equipartition, Owen et al. ( 2000) 1 find B-field 
strengths on the order of 7 - 10 //G, which are consistent with 
our lower limits. However, in cluster cores, where the particles 
and magnetic fields have different origins and evolutionary his- 
tories, the validity of the equipartition condition is not obvious. 
To measure a magnetic field of 10 yuG using X-ray observa- 
tions would require the ability to determine with significance 
a power-law flux 50 times lower than our upper limit, which, 
within the current calibration accuracies of ~10% is still out of 
reach. Using the Faraday rotation of the halo radio source in 
M87 _Dennison_ (1980) finds a magnetic field in the halo of 2.5 
juG. While the magnetic field deduced using the Faraday rota- 
tion measure is the average along the line of sight of the product 
of the magnetic field and the gas density, the lower limit of the 
magnetic field determined here depends on the volume aver- 
age of the relativistic electron density and on the square of the 
magnetic field over the emitting region. The fields determined 
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using these two methods can be quite different, and as shown by 
Goldshmidt & Rephaeli ( |1993| l the magnetic field determined 
using the radio and IC X-ray flux densities is in general smaller 
than the value determined using the Faraday rotation measure. 
Sanders et al. ( 2005 | l found a non-thermal X-ray emission 



component in the core of the Perseus cluster extending to a 
radius of ~75 kpc. Assuming this emission to be due to IC 
scattering of the cosmic microwave background and infrared 
emission from NGC 1275, they mapped the magnetic field in 
the core of the cluster. Within the radius of ~10 kpc they found 
a magnetic field between 0.5-3 fiG, while at larger radii their 
inferred magnetic field decreased to a value of ~0.1;uG. The 
higher values in the core are mostly consistent with the lower 
limits found in the core of M87, but the infeiTed magnetic field 
at larger radii is significantly smaller than the lower limits we 
found at similar radii in M87. 

For the radio relic in Coma, combining the radio data with 
XMM-Newton observations, [Feretti & Neumann |p006[ l found 
a lower limit on the magnetic field of 1.05 fiG, which is con- 
sistent with our lower limits found for the radio arms. The de- 
duced magnetic field in Coma based on the detection of tails 



of hard X-ray emission by BeppoS AX ( Fusco-Femiano et al. 
2004) 1 is ~0.2 fiG. However, the cluster region probed by the 
large field of view of the BeppoSAX PDS is much larger than 
the core of Virgo investigated in this paper. 

8. Conclusions 

We used deep XMM-Newton observations of the M87 halo to 
characterize the spatially resolved temperature structure and 
the chemical composition of the multi-temperature gas asso- 
ciated with the inner radio lobes. We found that: 

- Compared to a simple two-temperature fit, we obtain a bet- 
ter and more physical description of the spectra using a 
model which involves a continuous range of temperatures 
in each spatial bin. The range of temperatures of the mul- 
tiphase gas spans between ~0.6-3.2 keV. The cooling time 
of this gas is longer than the age of the lobes, indicating that 
no significant amounts of gas cooled out from the temper- 
ature domain observable in X-rays. If the multi-phase gas 
is distributed in many small spherical blobs with typical 
blob sizes of 0.1 kpc, similar to the widths of the filaments 
seen in the Chandra images ( Forman et al.|2007) , the ther- 
mal conduction must be suppressed by a factor of ~30-100. 
Such suppression may be provided by the magnetic fields. 

- We find a correlation between the amount of gas cooler than 
1 .5 keV and the metallicity, indicating that the cool gas is 
more metal-rich than the ambient halo. Extrapolating the 
linear fit between the cool gas fraction and Fe abundance, 
we estimate the Fe abundance of the cool gas to ~2.2 so- 
lar. This suggests the key role of the AGN in transporting 
heavy elements into the intracluster medium by uplifting 
cool, metal-rich gas from the galaxy. The abundance ra- 
tios of O/Si/S/Fe in and outside the X-ray arms are similar, 
indicating that the metals outside the multi-temperature re- 
gion, which were (according to our model) uplifted from 
the galaxy by previous AGN outbursts, must have been 



transported after the last major epoch of star formation, 
when contributions of SNll to the chemical budget of the 
galactic gas became negligible. 

- We estimate the mass of the gas with a temperature below 
1.5 keV (the smallest value in the deprojected profile of 
the hot M87 halo) to be about 5 x 10** Mq. Its chemical 
composition indicates that it originates from a mixture of 
ICM and stellar mass loss enriched with Type la supernova 
products. The amount of mass from stellar winds and SNela 
represents 15-50% of the total mass of uplifted gas. The 
time required to produce the observed metals in this gas 
is a;30-110 Myr, suggesting that the uplift of cool gas by 
AGN radio bubbles is a relatively rare event. 

- We put upper limits on possible non-thermal X-ray emis- 
sion from M87, and combining it with the 90 cm radio 
maps, we put lower limits of around 0.5-1.0 /iG on the mag- 
netic field strength. 
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